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DECLAR 
BAS 


; BASIC matrix multiply 


ATIONS 
TLMUL = Multiply 2 arrays giving a third 
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TITLE BASSMAT_MUL ; BASIC matrix multiply | 

SIDENT /1-021/ 3 File: BASMATMUL.MAR Edit: $8L1020 
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ae 
FACILITY: BASIC code support 


ABSTRACT: 


This module multiplies 2 arrays of any dtype and stores the result in a 
third array of any dtype. 


ENVIRONMENT: User Mode, AST Reentrant 


AUTHOR: R. Will, CREATION DATE: 11-Jul-79 
MODIFIED BY: 
++ 


SS SS SS SSS SSS SSSSSSSSssesssss 
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1-001 = Original 
1-308 - Change MTHSDFLOOR_R1 to MTHSDFLOOR_RS. JBS ab peat hy 
1-005 - Add check for Illegal Operation error. RW 28-S pt 
1-004 - Set IV bit in mask to signal integer overflow. RW 2-Oct-79 
1-005 - Redo scal HE Ri -Dec- 
0 1-006 - Change MTHS FLOOR_RS to MTHSDINT_R4. JBS 19-DEC-1979 
0 1-007 - Fix test for ‘same array’ for virtual. RW 15-Feb-1980 
3 1- 38 - Add support for byte, g and h floating. PLL 17-Sep-81 
1- - More modifications for new data types. PLL 24-Sep-81 
0 1- 19 - Changed shared external reference to G* RNH 25-Sep-81 
0 1-011 - Substitute a macro for the cette to the array fetch and store 
0 routines. This should speed th ngs up: PLL 9-Nov-81 
; 1-012 - gorres’ rya~t tee expression in the FETCH and STORE macros. 
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; BASIC matrix multiply 


: : 
8 «8 
000 64: 
3 
0000 $9 : 
0000 68: 1-021 - 
0000 69 ; 


BSE=198e 10:80:28 eupmewemaine mar: 1 


~ not store an hfloat Lonent , ” *. 15-Feb-82 
As 


t List macro expansions. 
pont pr : 


oblem 
cbenge _— 
code 


Change own storage to stac 


h stack t¢ nolned” by oane 0 PLL Fm oleh» 
ro instruct! ons at STORE frLoat: vee nar et 2 
to support errays of descriptors. 1982. 
storage. LEB 9-J jul 1 


Allow gfloat fosusts to HH} groves in a double at and 


vice versa. 
fix einer t 
MOL 15-Oct- 
Use G* for 


PLL 7-Oct-19 
types in ot word*tlong, and anything*thfloat. 


ALL externals. SBL 16-Nov-1982 
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Rete BECLARATIONS te 66-SEP-1984 13:40:58 YBASRTL. SREIGASMATHUL..MAR: 1 . 3) 
f ~SBTTL DECLARATIONS 
: ; INCLUDE FILES: 
fg 
7 SOSCDEF 3; define descriptor offsets 
4 SSF DEF 3; use to get scale 
0 8 3 
0 > EXTERNAL DECLARATIONS: 
33 ; -DSABL GBL Prevent undeclared 
8 3 4 symbols from being 
00 5 automatically 9'g al. 
0000 § -EXTRN BASSK_ARGDONMAT signalled if al blocks 
443 not present in array desc 
000 8 or dimct = 0 
0000 34 -EXTRN BASSK_DATTYPERR signalled if dtype of arra 
0000 0 isn't word Long float double 
0000 91 -EXTRN BASSK_MATDIMERR ; signalled if # of dims on any 
0000 38 error sn‘t 
0000 9 -EXTRN BASSK_ARRMUSSAM ; signalled if cols of srci not 
0000 94 = rows of src 
0000 -EXTRN BASSK_ILLOPE ; signalled if dest matrix is 
0000 same as either src matrix 


95 
96 

0000 4 -EXTRN BASSSTO_F 
99 
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A-B_R8 ; array element store for byte 
0000 eEXTRN BASSSTO_FA_W_R ; array element store for word 
0000 -EXTRN BASSSTO_FA_L_R ; array element store for long 
0000 100 ~EXTRN BASSSTO_FA_F_R ; array element store - float 
0000 1) -EXTRN BASSSTO_FA_D_R array element store - double 
0000 1 ¢ -EXTRN BASSSTO_FA_G_R8 array element store - gfloat 
0000 13 eEXTRN BASSSTO_FA_H_R8 array element store - hfloat 
000 104 eEXTRN BASSFET_FA_B_R array element fetch - byte 
000 105 -EXTRN BASSFET_FA_W_R array element fetch - word 
00 196 -EXTRN BASSFET_FA_L_R array element fetch = long 
4 10 ~EXTRN BASSFET_FA_F_R array element fetch - float 
00 108 -EXTRN BASSFET_FA_D_R array element fetch - double 
0000 109 ~-EXTRN BASSFET_FA_G_R array element fetch = gfloat 
000 110 -EXTRN BASSFET_FA_H_RB array element fetch = hfloat 
One 111 -EXTRN BASSMAT_REDIA check if redimensioning of 
00 11 dest array is necessary, if 
000 11 so, do it 
900 114 -EXTRN BASSSSCALE_R1 scale for double procision 
115 ~EXTRN MTHSDINT_R& routine to integerize double 
44 -EXTRN BASSSSTOP signal fatal errors 
0 11 -EXTRN BASSFETCH_BFA 
00 118 -EXTRN BASSSTORE—BFA 
0 11 
0 120 ; 
. 121 ; MACROS 
000 438° 
0 124 ; SBASSMAT_MUL multiply loop algorithm, see next page 
125 ; FETCH fetch an element from an array 
: § 3 STORE store an element into an array 
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00000000 
00000004 
00000008 
0000000C 
00000010 
00000014 
00000018 
0000001C 
0000002C 
00000034 


0000004 
0000048 
000004 
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So 


: PSECT 
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. 


EQUATED SYMBOLS: 


upper_bound_k = 0 


Lower _bound_k = 
upper_bound_j = 
lower _bound_j = 
upper_bound_i = 
current_j = 20 
current_i = 24 
current_sum = 28 
scale = 44 

srcl = 52 
value_desc ri 66 


wwronoen 
OoOontwos 


> OWN STORAGE: 


DECLARATIONS: 


12 
16 


-PSECT _BASSCODE PIC, USR, CON, REL 
EXE, RD 


, NOWRT. LONG 
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stack offset for semporary 
for upperbound of inner loop 
stack offset for sengerery 
for lowerbound for innerloop 
stack offset for songerory 
for upperbound of middle loop 
stack offset for temporary 
for lowerbound of middle loop 
stack offset for penperery 
for upperbound of outer loop 
stack offset for Sonperery for 
current value of middle loop 
stack offset for temporary for 
current value of outer loop 
stack offset for temporary for 
puaetag to get element 
otect offset for temporary for 
scale 
place to store element 1 while 
element 2 is fetched 


desc offset if sub 
desc offset if 1 sub 
desc offset if sub 
desc offset if sub 
desc offset if sub 
desc offset if sub 


SHR, = 


e 
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| BASRTL.SRCJBASMATMUL .MAR; 1 


z¢ 

; This macro contains the coping mechanism for accessing all elements of 

; an array. It also contains all the logic for all the combinations of data 

3; types and scaling. A macro is used to make it easy to maintain the parallel 
; code for all the different data types. 


7 -MACRO SBASSMAT_MUL srcl_dtype, src2_dtype ; multiply algorithm 


3¢ 
: Get scale so if any of the arrays is double, the scale will be there 


pass FP to get scale 

get scale in RO & R1 

call a BLISS routine because 
the frame offsets are only 
defined for BLISS 

store the scale 


MOVL SFSL_SAVE_FP(FP), RO 
JSB 


MOVD RO, scale(SP) 


3+ 

3; Loop eorours all the rows of the destination matrix. 

; Row and column upper and lower bounds have been initialized on the stack. 
; Current row (current_i) has been initialized to its lower bound. 


LOOP_1. ‘arc! dtype srce.dtypg y 
MOVL ower_bound_j(SP), current_j(SP) ; initialize current_j 
; to lower_bound of 


00 3+ 
00 ; Loop through all the elements (columns) of the current row of the destination 
00 ; matrix. Current column (current_j) has been initialized to its lower bound. 
80 ; Column upper bound is on the stack (upper_bound_j). 
00 : 
00 LOOP_J_"srcldtype’ src? _dtype 
00 MOVL ower _bound_k(SP), R11 : initialize current_k (R11) to 
Bp ; lower_bound of k 
0 CLRQ current_sum(SP) ; a CLRQ will set the temporary 
0 ; to 0 for all possible dtypes 
CLRQ current_sum+8(SP) ; sum could be hfloat (but don't 
- ; use h instruction here) 


Loop thr ugh. summing the products of each element of the ith row of srcl 
and the it column of srcé. current_i and current_j are on the stack. 
Source array element pointer (current_k) has been initialized in R11. 
Distinguish array by data type so that the correct fetch routine can 
retrieve the data, the correct multiply and add can be done and the correct 
store routine can be called. 


ad Be Se Ge Se Sete es 


OOP_K_'srcl_dtype'src2_dtype': 


g¢ 
3 Get the data from the first source array 
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MOVL srcl_matrix(AP), RO ; pointer to Ist src array 


ix multiply i oj tb 9 93:30:39 arth i Macro v04-00 Page (d) 


or? 


Pate a ee tae ee 
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MOVL C ‘vag tat R1 3 
MOVL R11, _R 

FETCH “grcl_gtype" 
MOV'srcl_dtype + srci(SP) 


+ 
; Get the data from the second source array 


MOVL srecmatrix(AP), RO 3 
MOVL R11, R1 : 
MOVL current j(SP), R2 5 
FETCH src2_dtype' : 


+ 
If the data types of the 2 source arrays is the 
Ife data type. 


f either o 


IF IDN srci_dtype, src2_dtype ; 

IF IDN srci_dtype, D 3 

MULD srci(SP), RO ; 

DIVD2 scale(SPs, RO : 

CMPD scale(SP), #1 3 

BEQL : 

JSB G*MTHSDINT_R4 3 

1$: ADDD2 =RO, current_sum(SP) : 


olFF 
MUL'srcl_dtype' srci(SP), RO 
0 tS ha al 


° a 
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1OD OD G9 G9 G99 GD G9 GD. O9 SII DDS DD TI EES BS EE EMA 


current row 


; current col 
; fetch data from srcl array 


store the Ist array element 


pointer to 2nd src array 
current row 

current col 

fetch data from src2 array 


same, do the arithmetic in 


Else convert the data to a common type and quitioly and add. 
the source elements is double, descale and multiply. Sc 
ntegerize the product before adding it to the sum. 


ale and 


; src arrays are 

; same data type 

; both sources are double 

; multiply 

; remove extra scale from result 
; is the scale 0? 

; yes. do not integerize 

; no, integerize 

; add to sum 

; sources same dtype, not double 


multiply the source elements 


RO, current_sum(SP) ; add product to current sum 


end of same dtype cote 
ype 


IF 3 src arrays different 
IF IDN srci_dtyp . 2 ; source 1 is hfloat 
CVT gree dtyee 0, RO 3; cvt array2 to hfloat 
MULH srci(SP), RO ; mult the elements 
mt RO, current_sum(SP) ; add product to current sum 
IF IDN src2_dtype, H ; source 2 is hfloat 
vt" rel dtype’H srci(SP), srcl(SP) ; cvt array! to hfloat 
RUL HG srci(SP), RO ; mult the elements 
anon RO, current_sum(SP) 3; add product to current sum 
a | IDN srci_dtype, G ; source 1 is gfloat 
0 IF IDN src2_dtype, D : special case if g & dbl 
1 Divd2 scale(sP), RO ; descale src 
¢ CVTDOH =RO, R 3; cvt src2 to hfloat 
CVIGH = srci(SP), srci(SP) i; cvt srcl to hfloat 
4 prety srci(SP), ; mult the elements 
5 ADDH RO, current_sum(SP) ; add product to current sum 
§ LFF 3; src2 any type but dbl 
CvT' ree_dtype'G RO, RO i cvt srce to gfloat 
8 RL GS srci(SP), RO ; mult the elements 
9 ADDG RO, current_sum(SP) ; add product to current sum 
90 -ENDC : end of srci gfloat 
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o LFF 
lf IDN srcé_dtype, G ; source 2 is gfloat 
If ION src sot yee Dd i special case gfloat & dbl 
DIVD2 scale(SP), srcl(SPS ; descale src 
CVTOH = srci(SP), srci(SP) ; cvt srcl to hfloat 
CVTGH RO, 3; cvt src2 to hfloat 
prey srci(SP), RO 3; mult the elements 
ADOH RO, current_sum(SP) 3; add product to current sum 
LFF ; src2 gfloat & srci not dbl 
CvtT' rel dtype srcl(SP), srcl(SP) ; cvt srcl to gfloat 
MULG srci(SP), RO ; mult the elements 
ADDG RO, current_sum(SP) ; add product to current sum 
at de ; end of src2 gfloat 
IF IDN srci_dtype, D ; source 1 is double 
; don't have to worry if src2 
; is gfloat because we already 
i: checked 
CVT'src2_dtype'D RO, RO ; cvt array2 to double & save 
; no scaling needed because in 
3 mult iely scale is in 2nd src 
MULD2 srci(SP), RO 3; multip Y 
JSB G*MTHSDINT_R4 ; integerize 
ADDD2 = RO, current_sum(SP) ; add product to current sum 
o LFF ; Ist_array not double 
° IDN src2_dtype, D : is 2nd sre double 
; yes, make srcl double & save 
CVT'srcl_dtype'D srci(SP), srcl(SP) ; make Ist src double 
; no scaling needed because for 
3 mult ety only 1 number is 
; needed to be scaled. 
MULD2 srci(SP), RO 3 compute the product 
JSB G*MTHSDINT_R4 ; integerize 
ADDD2 =RO, currenft_sum(SP) 3; add product to current sum 
LFF 3 no double operands try float 
ti IDN srcl_dtype. F : is Ist element float 
CVT gree. dtyee.t 0, RO ; make 2nd element float 
reed srci(SP), RO ; multiply the elements 
ADDF RO, current_sum(SP) ; add to current sum 
LFF 3; ist_array not float 
IF DN src2_dtype, F ; is 2nd array float 
CvT rel dtypert srci(SP), srci(SP) ; yesomake Ist element float 
RULES srci(SP), ; multiply the elements 
ADDF RO, current_sum(SP) ; add to current sum 
LFF 3; no double or float, try long 
* IDN srcl_dtype. L 3; is Ist array long 
CvT' ree_dtype rt 0, RO 3; make 2nd element long 
RULLS srci(SP), RO ; multiply elements 
Apot RO, current_sum(SP) ; add product to current sum 
IF IDN src2_dtype, L ; src2 is long 
CvT° rel dtype t srci(SP), srci(SP) ; cvt srcl to long 
RULE srci(SP), RO ; multiply the elements 
Appt RO, current_sum(SP) 3; add product to current sum 
ti IDN srcl_dtype. wd ; srecl is_word 
CVT'src2_dtype'w 0, RO 3; cvt src 


BASSMAT 
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BASIC mat 
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+ 


it is the 
Otherwise, store the sum in the destination arra 
(pass pointer to dest in R10 and pointer to stac 
and continue with next destination element. 


5$: 


ix multiply 
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pert | srci(SP), RO 3; mult the elements 

tii RO, current_sum(SP) ; add product to current sum 
aU ION srcl_dtype. 8 ; srcl is byte 

CvT' ree dtype 8 » RO 3 cvt src2 

MULB srci(SP), RO 3; mult the elements 

ADDB RO, current_sum(SP) ; add product to current sum 
LFF ; last case - src2 must be b 
CvT' rel dtype 8 srci(SP), srci(SP) ; cvt srct 

nULBS srci(SP), RO ; mult the elements 

a ttt: RO, current_sum(SP) 3; add product to current sum 
~ENDC 

eENDC 

wENDC 

~ENDC 

-ENDC 

-ENDC 

eENDC 

eENDC 

~ENDC 

eENDC 

-ENDC 

eENDC 


ast product of the sum. 


If 


Have multiplied next set of elements and added it to current sum. 


See if 


not continue with current sum 


INCL R11 

CMPL R11, upper_bound_k (SP) 

BGTR 5$ 

BRW LOOP_K_‘srcl_dtype'src2_dtype 
MOVL SP, RS 

IF IDN srci_dtype, src2_dtype 
Bipy DEST_CASE_‘srci_dtype' 

c IF IDN srci_dtype, H 

od 1 DEST_CASE_H 

vIF IDN src2_dtype, H 

4 DEST_CASE_H 

v1F 1DN srcl_dtype, G 

IF IDN src2_dtype, D 

BSBW DEST_CASE_H 

J IFF 

BSBW DEST_CASE_G 

-ENDC 


by calling a subrout ine 
in R5) 


; get nex 


: t K 
3; see if last product in sum 


no, continue inner loop 


finished inner loop so sto 
point to temps 
src arrays are 
same data type 
go cvrt to dest type and s 
src arrays different dtype 
source 1 is hfl 


source 2 is hfloat 


lLoat 
special case gfloat & dbl 
wps done in hfloat so cvt 
from hfloat to dest type 
gfloat & all other dtypes 


source 1 is gfloa 
d 


Page (d 


yte 


re 


tore 


oat 
cvt from hfloat to dest type 


cvt from hfloat to dest type 


Cvt from gfloat to dest type 
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405 -IFF 
4 $ lf IDN src2_dtype, G i; source 2 is gfloat 
4 If IDN srci_dtype, D i; special case dbl & gfloat 
00 4 BSBw DEST_CASE_H i ops done in hfloat so cvt 
00 4 i; from hfloat to dest type 
00 410 FF ; gfloat & all other dtypes 
00 411 BSBW DEST_CASE_G i; Cvt from gfloat to dest type 
00 216 ~ENDC 
B85 41 1FF 
00 414 IF IDN srci_dtype, D ; source 1 is double 
000 415 BSBW DEST_CASE_D 3; cnvrt from double to dest type 
§808 $18 ; and store 
00 4 i ; (note that we don't have to 
0000 418 ; worry about dbl & gfloat here 
0000 419 ; .because it was handled above) 
0000 420 o LFF i; Ist_array not double 
0000 421 IF IDN src2_dtype, D ; is 2nd src double 
0000 4 § 3; yes, make srcil double & save 
0000 4 BSBW DEST_CASE_D i; cnvet from double to dest type 
0000 424 ; and store 
0000 425 LFF ; no double operands try float 
0000 426 IF 1DN srci_dtype, F i is 1st element float 
0000 427 BSBw DEST_CASE_F 3; cnvrt from float to dest type 
0000 428 3 and store 
0000 429 oIFF i Ist_array not float 
0000 430 IF IDN src2_dtype, F 3 is 2nd ervey float 
0000 431 BSBW DEST_CASE_F i cnvrt from float to dest type 
464 $36 ; and store 
000 43 -IFF 
0000 434 IF IDN src* dtype, L ; source 1 is long 
0000 435 BSBW DEST_CASE_L 3; cnvrt from long to dest type 
0000 436 3; and store 
0000 437  IFF 
0000 438 IF IDN src2_dtype, L ; source 2 is long 
0000 439 BSBW DEST_CASE_L i cvt from long to dest type 
0000 440 FF 
0000 441 IF IDN srci_dtype, W ; source 1 is word 
0000 $66 BSBW DEST_CASE_W i; cvt from word to dest type 
0000 44 LFF 
0000 444 IF IDN src2_dtype, W ; source 2 is word 
0000 445 BSBW DEST_CASE_W i; cvt from word to dest type 
0000 446 3; byte and any other data type 
0000 447 ; would've been caught by one 
6008 723 ; of the above cases, and byte 
0 44 ; & byte is handles by case 1 
it 450 FF 
8 00 451 If IDN srci_dtype, % ; source 1 is byte 
000 $26 BSBW DEST_CASE_B 
000 45 IFF . 
000 454 BSBW DEST_CASE_B i only thing left 
000 455 -ENDC 
444 $28 -ENDC 
00 45 -ENDC 
00 458 ~ENDC 
Bh 459 ~ENDC 
0 460 -ENDC 
0000 461 ~ENDC 


ix multipl 15-SEP-1984 23:47:5 AX/VMS Macro v04-00 Pa 
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46 -ENDC 

re eENDC 

464 -ENDC 

465 -ENDC 

46 -ENDC 

46 -ENDC 

re -ENDC 

3 

470 ;+ 

471 ; Have stored that element. Now see if it was the last column. If not, 
ti 3 continue with the next column. Otherwise continue to next row. 

47% ° 

475 INCL current_j (SP) 3; get next column 

$78 asa suerent.. (sp), upper_bound_j(SP) ; see if last column done 
$78 BRW LOOP_J_*srcl_dtype'src2_dtype* ; no, continue inner loop 
480 ;+ 

2B) ; Have completed entire row. See if it was the last row. If not, 

rt} 3 continue with next row. 

oF 

485 208: INCL current_i(SP) 3; get next row 

+38 veal Fucrent 15h? upper_bound_i(SP) ; see if last row done 
r+) BRW LOOP_I_"srcl_dtype'src2_dtype’ ; no, continue outer loop 
rs 108: RET : yes, finished 


BASSMAT_MUL ; BASIC matrix multipl 15-SEP-1984 23:47:5 AX/VMS Macro v04-00 Pa 
1-021 BASSMAT_MUL - Multiply 2 arrays giving g-$ep- 1986 93:40:35 EBASRTL. SRCIBASMATMUL..MAR; 1 a 
| $38 bce ~SBTTL BASSMAT_MUL = Multiply 2 arrays giving a third 
496 ; FUNCTIONAL DESCRIPTION: 
00 498 ; Multiply 2 arrays giving a third. Signal an error if the er and 
| 0 499 ; lower’ bounds (exe lading 0) for columns in srci_matrix do s not equal 
00 0; the upper and lower bounds (excluding 0) for rows in src2_matrix. 
00 13 An error will also be signalled if any of the three matrices does not 
000 ; 3 have a DIMCT of 2, or if DSCSA_POINTER in either srcl_matrix or 
Boh 3 srcé_matrix is the same as DSCSA_POINTER of dest_matrix. 
00 4; Redimension the output to have a lower bound of 0 for both dimensions, 
4 8? 3 and an upper bound for rows equal to the upper bound for rows for 
0 $ ; srci_matrix, and an_upper bound for columns equal to the upper bound 
000 07 ; for columns for src2_matrix. Initialize all the necessary 
0000 205 3 looping information on the stack. Conversions may have to be done 
444 09 ; so that the sources are the same data type, so divide 
000 239 ; the looping portion according to the data nypes. Conversion to the 
0000 11 ; correct destination data type will be done by a JSB to a routine, 
444 21§ ; instead of multiplying the number of possible combinations by 4. 
909 314 ; CALLING SEQUENCE: 
44 316 ; CALL BASSMAT_MUL (srcl_array.rx.da, src2_array.rw.da, dest_matrix.wx.da) 
9009 18 : INPUT PARAMETERS: 
00000006 0000 520 ° srcl_matrix = 4 
00000008 0000 521 src2_matrix = 8 
0000 5 § $ 
B00 2 ? : IMPLICIT INPUTS: 
44 5 2 ; Scale from the callers frame to scale double precision. 
9000 ; f : OUTPUT PARAMETERS: 
0000000c 0000 9° dest_matrix = 12 
0000 0; 
000 1 ; IMPLICIT OUTPUTS: 
000 ¢ 3 
ef Sa 
0000 2 5 > FUNCTION VALUE: 
44 § ; COMPLETION CODES: 
$600 3 8 ; NONE 
44 9; 
ts rt ; SIDE EFFECTS: 
0000 $5 : This routine calls the redimensioning routine and the array element 
000 43 ; fetch and store routines and therefore may signal any of their errors. 
000 44; It may also signal any of the errors Listed in the externals section. 
000 45 ; It may also cause the destination array to have different dimensions. 
000 84? ; 
000 $8 : 
GFFC 000 4 ENTRY BASSMAT_MUL, “M<R2,R3,R4,R5,R6,R7,RB,R9,R10,R11,1V> 
0002 50 
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"B-SEp=198e 10:40:38 EBASRTC. SRESBASMATAOL mar: 1 


REGISTER USAGE 
RO -R cootroves by store routines 
e 


R net us 

R10 pointer to dest matrix descriptor sexcept for double in which 
case R10 is part of double value R9-R1 

R11 current value of inner subscript 


oa 


Put routine arguments into registers for ease of use. 
If block 2 of array descriptor (multipliers) is not present then error. 


fall into error 


MOVL srcl_matrix(AP), R2 ; ptr to srcl array descr 
BBC #DSCSV_FL_BOUNDS, DSC$B_AFLAGS(R2), ERR_ARGDONMAT 

z exit Tf block 3 not 

: po in descriptor 
CMPB DSCSB_DIMCT(R2), #2 : 2 dimensional? 
BNEQU ERR_MATDIMERR :; if not, error 
MOVL src2_matrix(AP), R3 ; ptr to src2 array descr 
BBC #DSCSV_FL_BOUNDS, DSCSB_AFLAGS(RS), ERR_ARGDONMAT 

; exit Tf block 3 not 

; present in descriptor 
MOVL dest_matrix(AP), R10 3 pointer to dest descriptor 
CLRQ -(SPJ 3 Save space for VALUE_DESC 
CLROQ -(SP) 3 AND 
CLRQ -(SP) 3; DATA 
CLRD -(SP) ; space for srcl element 
CLRD -(SP) ; it may be hfloat 
TSTD (SP) 3; save space for scale 
CLRO -(SP) ; reserve space to save sum 
CLROD (SP) ; _possible hfloat sum 
CMPB DSCSB_DIMCT(R3), #2 5 dimensional? 
BEQLU INIT TwO SUBS : if 2-D continue, else 


ERR_MATDIMERR: 
PUSHL 


#BASSK_MATDIMERR Signal error, src arrays 


CALLS #1, G*BASSS$STOP : don't have same # dimensns 
ERR_ARGDONMAT : 
PUSHL #BAS$K_ARGDONMAT ; signal error, 0 for dimct 
CALLS #1, G*BASS$S$STOP : or block 2 or 3 absent 
ERR_ARRMUSSAM: 
PUSHL #BASSK_ARRMUSSAM ; Signal error, src arrays 
CALLS #1, G*BASS$S$STOP 3 same bounds 
ERR_ILLOPE: 
PUSHL #BASS$K_ILLOPE : Signal error, dest array is 
CALLS #1, G*BASSS$STOP ; aS one of source arrays 


Zé 
; There are 2 subscripts. Check and redimension the destination array if 
3 necessary. Put the upser bound for both subscripts on the 


; stack and make sure that the lower bound for both subscripts will start 
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20 A328 A2 CMPB Ss dsc $l_u2_2(R2), dsc$l_ul_2(R3) dges src array have the same 
nd upper bound as src2 array 
has for — upper bound 
no, erro 
does sre] array have the same 
2nd Lower bound as src2 array 
bound 


has for Ist lower 


"OCCT @OMwow 


BNEQU ERR {Rens 
CMPB dsc $L_l2_2(R2), dsc$i_l1_2(R3) 
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no, error 
qne upper bound 


BNEQU ERR trueagts 
st upper bound 


dest array pointer 


CALLS G*BASSMAT redimension destination 


betes CLASSCRTON. MIDSCSK_CLASS_A : is array virtual? 
10 BNEQU TUAL _ SAME ; yes go check virtual 
3 dest not virtual even if 


: ere is virtual pointer will 
as ° and won't match 
$A reserenthe) DSCSA _POINTER(RI . - is dest same as srcl 
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10 A324 A2 
D 
a 
OA 
00000000‘ GF 8 
04 
04 AA OG A D1 CMPL DSC 
C 13. 0090 BEQLU ERR_ICLO yes rror 
04 AA 04 A3 D1 009 CMPL DSCSA POINTER(RS), DSCSA _POINTER(R 0)° 3 “te dest same as src2 
C 13 9 BEQLU ERR_ICLO 3 yes. error 
2 11 0099 6 BRB INIT_ Stack 
98 o3 VIRTUAL_ SAME: 
04 03 A2 94 9B 6 CMPB DSCSB_CLASS(R2), #DSCSK _CLASS_A ; ~ srcl virtual? 
0 13 OO9F 63 BEQLU 1$ cant be same 
FC AA FC ») Boag o37 Hag a LOGUNIT(R2), DSCSL stocunrinfBs | is dest same as src1? 
FB AA FB A2 O01 OOAS $38 CMPL DSCSL_BYTEOFF(R2), DSCSL _BYTEOFFiR10) ze is dest same as src1? 
BRA8 640 (ch . logunit and byteoff) 
A 13 AD 641 BEQL ERR_ILLOPE a rror 
04 O3 A 91 OAF 64g 18: CMPB pecsee CLASS(RS), #DSCSK_CLASS_A - sree, virtual? 
OE 13 008 64 BEQLU INI n't be same 
L : uN L_LOGUN ; 1s dest same as srcl? 
FC AA FC 4 0 #4 ote Ha bscsC LO og IT(R2), DSCS$L_LOG rings $ 4 1? 
FB AA FB AS D1 goer $68 CMPL pscsc BYTEOFF.(R3), DSCSL -BYTEOFF (R10) : is dest same as src2? 
98 13 pct o et state ERR_ICLOPE 3 yes error 
1C ag DD bes og PUSHL dsc$l_l1_2(R2) ; initialize current_i counter 
0 14 org 650 BGTR 1$ ; not row 0 or neg 
6€ ot dO C $2) MOVL #i, (SP) 3 start with 1 
E D4 O00CB 6 é 1$ CLRL 3s = ¢ §p) : Save space for current j 
OA DD cD «65 PUSHL dsc$t vu} 3 (Re) 3 upper_bound_ 
4A DD 00D 654 PUSHL sc$l_l2_2(R3) ; lower_ =boundj 
0 14 D 655 GTR $ 3 not row 0 or ony: do cols 
6— 01 00 D $26 VL #1, (SP) 3 start with r 
8 AS DD D8 657 2$ PUSHL qsc$t ug 5 (3) ; upper_bound_ 
SA DD poe $28 PUSHL gecSl_t _2(R2) 3 lower— spound 
0 14 D 65 3 not col 0 oF neg 
6— 01 00 E 660 MOVL #1, (SP) 3 start with k=1 
8 A2 DD E 661 3$: PUSHL dsc$l _u2_2(R2) 3 upper_bound_k 
Fe ee8 ss 
Oe 664 ; Algorithm now differs according to data types 
| 
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; BASIC matrix multipl 15-SEP-1984 23:47:5 AX/VMS Macro v04-00 Pa 14 
BASSAAT RUC’ = RULTIDLY 2 arrays giving ‘SosePolobe 1O;90;2% EBASHTe eae TASHAToR mar; 2% 18, 
E 5 ;- 
E $6 
+ SEPARATE _DTYPES: 
05 06 02 A2_ 8F O0E6 089 5$: CASEB DSCS$B DTYPE(R2), #DSCSK_DTYPE_B, #<DSCSK_DTYPE_D = DSC$K_DTYPE_B> 
7* QOEB 2 9 2s: -WORD BYTE- ; code for byte dtype 
F61° OOED 7 «WORD WORD-2$ ; code for word dtype 
3 * QOEF ore «WORD LONG-2$ ; code for long dtype 
* OOF 67 eWORD ERR_DATTYPERR-2$ ; quad not supported 
DB3" OOF 674 -WORD FLOAT-2$ z code for float dtype 
cpo9° : erg «WORD DOUBLE-2$ 3; code for double dtype 
‘7 77 ;¢ 
ore ers ; G and H floating fall outside the range of the CASEB. 
fF? 680° 
18 02 a2 91 OOF? 681 CMPB SCSB_DTYPE(R2), #DSCSK_DTYPE_G 
0 1 OFB 6 § BNEQ 
4C1 3 A ? Z BRw GFLOAT ; code for gfloat dtype 
1¢ 02 a2) 91 0100 685 38: CMPB «sd DSCSB_OTYPE(R2), #DSCSK_DTYPE_H 
0 \¢ 106 6 § BNEQ 4 
SB6F 3 i ? 5 BRW HFLOAT ; code for hfloat dtype 
18 O2A2 91 +4 689 4$ CMPB DSCSB_DTYPE(R2), #DSCSK_DTYPE_DSC 
08 is 100 690 BNEQ ERR_DATTYPERR 
52. O4 A D 10— 691 MOVL 4(R2), R2 ; R2 <== addr of descriptor 
D1 so Wa $36 BRB 5$ 3; CASE again on dtype in desc 
011 69 
0115 694 ERR_DATTYPERR: 
00000000'8F DD 0115 695 PUSHL #BAS$K_DATTYPERR : Signal error, unsupported 
00000000'GF 01 FB 0118 696 CALLS #1, G*BASS$SSTOP 3 dtype in array desc 
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15-SEP=1986 23:47:50 VAX/VMS Macro v06-00 Page 1 
5-$Eb=19be 10:90:24 EBASKTe eRe TbAMMATOCL.man:1 2% 1, 


» MDSCSK_DTYPE_B, #<DSCSK_DTYPE_D = DSCSK_DTYPE_8> 
; code for byte dtype 

code for word dtype 

code for long dtype 

quad not supported 

code for float dtype 

; code for double dtype 


18 4O2A 


SSN Oe 


CMPB DSCSB_DTYPE(R3), #DSCSK_DTYPE_G 
BRW § BYTE_TO_GFLOAT 

DSCSB_DTYPE(RS), #ADSCSK_DTYPE_H 
BRW = BYTE_TO_HFLOAT 

DSCSB_DTYPE(R3), #DSCSK_DTYPE_DSC 


MOVL 4(R3), R3 3 R3 <-- addr of descriptor 
bee BYTE 


BRW ERR_DATTYPERR 


; Now type of sourcel and source2 arrays are known. Use the macro to 
; generate the code for each case 
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; CASE again on dtype in desc 
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§ 6 and H floating fall outside the range of the CASEB. 
3 
; 
5 
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70 738 BYTE.TO_WORD: — SBASSMAT_MUL 
BF 739 
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8780 


74] BYTE_10_LONG: 


15S 


SBASSMAT_MUL 


pai9be 10:80:28 
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BASSMAT_MUL BASIC matrix multiply 15-SEP-1984 $3797 +30 AX/VMS Macro v04-00 Page 19 
1-021 BASSMAT T MUL = Multiply 2 arrays giving 6=-SEP-1984 10:30:2 BASRTL.SRCIJBASMATMUL .MAR; 1 (5), 


0780 744 BYTE_TO_FLOAT: $BASSMAT_MUL 8, F 
0901 745 


L 10 | 
baSsmat mo mo mutt ipty 2 arrays giving "3 Pa} 984 19; $5; 38 YOASKIC. BAM SRATAUL.M MAR; 1 hee a) 


Q9D1 747 BYTE_TO_DOUBLE: SBASSMAT MUL B,D 
F 48 


41 
RASS TAT fut baSsha ae muttip BLY 2 arrays Pn _ Ep =}382 9; 45; 39 cre a tat BA PATROL « oMAR; 1 neat 2) 


BES £30 BYTE_TO_GFLOAT: SBASSMAT_MUL B, G 


RASSRAT UL ; BASIC matrix multiply 19-360" 1 988 $3347:30 AX/VMS Macro V04-00 Page 22 
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0E22 753 BYTE_TO_HFLOAT: SBASSMAT_MUL  B, H 
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BASSMAT_MUL ; BASIC matrix multipl 15-SEP-1984 23:47:5 AX/VMS Macro V04-00 Pa 23 
7-821 BASSMAT_MUL - Multiply 2 arrays giving 3-3 eo -1 ope $3: 40:35 EBASRTL. SREIBASMATMUL..MAR: 1 ” (5) 
104C 755 ;+4 
| : rt: i ; Sourcel array is a word array. Now differentiate on the source2 type. 
106¢ «6758 * 
05 06 02 a3 er 104C 759 WORD: CASEB DSCSB_DTYPE(R3), #DSCSK_DTYPE_B, #<DSCS$K_DTYPE_D = DSCS$K_DTYPE_B> 
02D" 1 ? 760 1$: «WORD WORD _TO_BYTE-1$ ; code for byte dtype 
£eh. 1 763 «WORD WORD _TO_WORD-1$ ; code for word dtype 
468" 105 r6¢ eWORD WORD TO LONG-1$ ; code for long dtype 
FOC4 105 76 -WORD ERR_BATTYPERR-1$ 3 quad not supported 
0689° 1928 764 «WORD WORD_TO_FLOAT-1$ ; code for float dtype 
OB8AA' : 28 163 «WORD WORD _TO_DOUBLE-1$ 3; code for double dtype 
198 eS $* 
105D 768 ; G and H floating fall outside the range cf the CASEB. 
105D 769 ;- 
1920 770 
18 02 ay 91 105D 771 CMPB DSCSB_DTYPE(R3), #DSCSK_DTYPE_G 
0 12 1061 ie BNEQ $ 
OaBC Ué31 Se? ot? BRw WORD_TO_GFLOAT ; code for gfloat dtype 
10 02 A391 1066 775 28: CMPB SCSB_DTYPE(R3), #DSCS$K_DTYPE_H 
03 12 106A £76 BNEQ 
OcDD ss é31 #3 ahs BRW WORD_TO_HFLOAT : code for hfloat dtype 
18 02 A3 91 106F 779 38 CMPB sd DSC SB_DTYPE(R3), #DSCSK_DTYPE_DSC 
06 12 1073 780 BNEQ ; 
53 04 A3 DO 1075 781 MOVL 4(R3), R3 ; R3 <-- addr of descriptor 
Dis + res BRB WORD 3; CASE again on dtype in desc 
F097, -31«107BSs 784 4S: BRW ERR_DATTYPERR 
107E 785 
107E 786 ;+ 
107E 787 ; Now type of source! and source2 arrays are known. Use the macro to 
1G7E 788 ; generate the code for each case 
107E 789 ;:- 
107E 790 
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ply 15-SEP-1984 ee 
ba SSMAT_MUL = Multiply 2 arrays giving 6=-SEP-1984 10:30 BASRTL.SRCJBASMATMUL .MAR; 1 (5) 
1 138 WORD_TO_BYTE: $BASSMAT_MUL W, B 
1290 


8 wots 


7? 
0 


{0:3 


P=1984 
b-seP= 1986 


155 
SBASSMAT_MUL 


795 WORD_TO_WORD: 
796 


89 


wT 
c- 


E11 


matrix multiply 15-SEP-19 


MUL 


= Multiply 2 arrays giving 6-SEP-19 
7 798 WORD_TO_LONG:  $BASSMAT_MUL Ww, 


Be {0:30 


L 


:23 


AX/VMS Macro v04-00 
BASRTL.SRCIBASMATMUL .MAR; 1 


a 


F 11 
BASSMAT_MUL ; BASIC matrix multipl 15-SEP-19 247:5 AX/VMS Macro v04-00 P 27 
1-021 BASSMAT_MUL - Multiply 2 arrays giving onsEP-198e $3:$0:35 UBASRTL. SRCJOASMATMUL MAR: 1 ses (5) 
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; BASIC matrix multipl 15S AX/VMS v04-00 Pa 31 
BASSMA MUL = Multiply 2 arrays giving 6-S pa} 984 Te 4: 38 YeASRiL.s L.S SSASMATMUL MAR; 1 we (5) 
: ;* ‘ 
: _Sourcel array is a longword array. Now differentiate on the source2 type 
7 "ae 
02 a3 er 7 ieee CASEB ty pr veg (Ro). #DSCSK_DTYPE_B, _M<DSCSK _DTYPE_D = pe _DTYPE_8> 
8 D® 1F7 : «WORD G_TO_BYTE-1$ code for byté dty 
4C* 17 «WORD LONG: O-BORD=1$ ; code for word dtype 
1 7 «WORD LON =1 ; code for long dtype 


; quad not supported 
; code for float dtype 
; code for double dtype 
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G and H floating fall outside the range of the CASEB. 
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02 a3 91 DSCSB_DTYPE(R3), #DSCSK_DTYPE_G 
03 12 8 BNEQ 
OABA 31 : BRW LONG_TO_GFLOAT ; code for gfloat dtype 
02 a3. 91 9 2s CMPB SCSB_DTYPE(R3), #DSCSK_DTYPE_H 
03 12 9 BNEQ 
ocobB sé 4 BRW LONG_TO_HFLOAT : code for hfloat dtype 
02 a3 91 9 837 3$ CMPB DSCSB_DTYPE(R3), #DSCSK_DTYPE_DSC 
0 12 9 8 BNEQ 
O46 A b0 9 8 MOVL 4(R3), R3 ; R3 <== addr of descriptor 
Di so > 8 BRB LONG ; CASE again on dtype in desc 
—16D 31 : : 4$: BRW ERR_DATTYPERR 
3+ 
A 8 : Now type of sourcel and source2 arrays are known. Use the macro to 
. ; ; _generate the code for each case 
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RASS TAT ut iaSSaat mo Mult 3 4 2 arrays giving "3 3Ep= 138s 19: 45: 33 YBASRTL. § BAS SRATRUL. MAR; 1 rte za 


gic? ees LONG_TO_WORD: S$BASSMAT_MUL L, W 


41 
baSSmat mo mor mut tie He 2 arrays whee. y* p= 1982 tf 45: 39 YBASKIL Bae SBASMA TOOL. MAR; 1 Aad i) 
$768 $36 LONG_TO_LONG: S$BASSMAT_MUL L.L 


N11 
BASSMAT_MUL BASIC matrix multiply 15-SEP-1984 $3! 247: +30 VAX/VMS Macro v04-00 Page 35 
1-021 BASSHAT MUL - Multip ply 2 arrays giving 6-SEP-1984 10:30:2 CBASRTL.SRCJBASMATMUL .MAR; 1 (5) 


2602 857 LONG_TO.FLOAT: SBASSMAT.MUL —L, F 


Pre en eee ne TEE @ Wigan 


ASSMAT_MUL 4,243! nm multiply 15-SEP-1984 7:5 AX/VMS Macro Vv04-00 Page 
; . - Multip ty 2 arrays giving goseP-1 984 9: $0: 38 EBASRIL. St SRE JBASMATMUL .MAR; 1 2, 


2823 560 LONG_TO_DOUBLE: SBASSMAT_MUL L. D 


sy 


cc 


BASSmA 


BAS! 


matrix multiply 
MUL 


- Multiply 2 arrays giving 


15-SEP-1984 
6-SEP-19 


ASA 863 LONG_TO_GFLOAT: SBASSMAT » 
tse $8? agi a , 


86 §0: 


47:5 
:30:2 


9 ort ond Macro v04-00 
BASRTL SRE IGASMATRUL. MAR; 1 


Pa 37 
a (5) 


bas 


BAS! 


SMA 


55 


C mat 
TLMUL 
C74 
ESE 


h) 
trix multiply 
= Multiply 2 arrays giving 


15-SEP-19 
6-SEP=19 


Be 79:50:25 


B66 LONG_TO_HFLOAT: SBASSMAT_MUL L, H 


AX/VMS Macr 
BASRTL.SRCJ 


ov 
BAS 


v04 
MA 


-00 
TMUL .MAR; 1 


w 8, 


—————— —— 
E 12 


BA SMAT_MUL ; BASIC matrix multipl 15-SEP-1984 23:47: AX/VMS Macro V04-00 Pa 39 
1-021 BASSMAT_MUL = Multiply 2 arrays giving g-SEp-} 984 $3:40:38 EBASRTL. SREIGASMATRUL MAR: 1 aie.” 
ESE 3¢ 
343 2 : Sourcel array is a floating array. Now differentiate on the source2 type 
E9E 872° 
0S 06 02 A3  8F £9 = FLOAT: CASEB DSCSB_DTYPE(R3), MDSCSK_DTYPE_B, #<DSCSK_DTYPE_D = DSC$K_DTYPE_B> 
9 D’ OEA 74 1$: .WORD FLOAT-TO_BYTE-1$ 3 code for byte dtype 
C* ZEA 75 «WORD FLOAT_TO_WORD-1$ ; code for word dtype 
398" EA 6 «WORD FLOAT_TO _LONG-1$ : code for long dtype 
D2 EA9 7 -WORD ERR_DATTYPERR-1$ 3 quad not supported 
besa EAB 8 .WO FLOAT_TO_FLOAT-1$ ; code for float dtype 
6' 4 a0 -WORD FLOAT_TO_DOUBL-1$ ; code for double dtype 
EAF 1 3¢ 
a Ho) ; G and H floating fall outside the range of the CASEB. 
EAF 884° 
91 2EAF 885 CMPB =: OSCS$B_DTYPE(R3), #DSCSK_DTYPE_G 
12 2€83 Hy BNEQ es 
31 EBS 887 BRW FLOAT_TO_GFLOA 
91 2688 889 28: CMPB =: OSC$B_DTYPE(R3), #DSCS$K_DTYPE_H 
12 2—€BC 890 BNEQ 3$ 
31 +41 +h BRW FLOAT_TO_HFLOA 
18 02 A3 91° 2EC1 398 3$: CMPB  =—s-: OSC $B_DTYPE(R3), #DSCSK_DTYPE_DSC 
06 12 2ECS 894 BNEQ = 4$ , 
53 04 A353 00 sete 895 MOVL 4(R3), R3 ; R3 <== addr of descriptor 
rn | cc. 538 BRB FLOAT ; CASE again on dtype in desc 
D245 31 ECD 43 4$: BRwW ERR_DATTYPERR 
EDO 900 : Now type of sourcel and source2 arrays are known. Use the macro to 
4 901 ; generate the code for each case 


1B 02 AS 
03 

OAB8 

1¢ ©6602 A3 
03 

OCD9 


--- 


BASSMAT_MUL 
1-021 


BASSMaT 


BASIC matrix multiply 
TLMUL = Multiply 2 arrays giving 


OEF 


EDO 904 FLOAT_TO_BYTE: $BASSMAT_MUL 


15-SEP-1984 
b-$6P=19 


84 23: 
84 10: 


47: 
30: 


:23 


VAX/VMS Macro v04-00 
CBASRTL.SRCJBASMATMUL .MAR; 1 


Page 


a“ 
Wo 


) 


er 


G1 
8 1 eos baSSRA wm muteipt 2 arrays ae "BSE 21386 19: 46: 3 YBASRTL. sacie JBAS SPATRUL.M MAR; 1 ne (> 


HH 307 FLOAT_TO_WORD: SBASSMAT_MUL F, W 


baSsmat mon’ mutt pty 2 arrays ne y = HH -138¢ $3: 47: 39 YAX/VMS- SRESBASRATA et Page 42, 
3298 310 FLOAT_TO_LONG: $BASSMAT_MUL fF, L 


adie 7 RTS Sa 112 


BASSMAT_MUL ; BASIC matrix multiply 15-SEP-1984 $3:47:39 AX/VMS Macro v04-00 Page 43 
1-021 BASSMAT_MUL = Multiply 2 arrays giving 6-SEP-1984 10:30:2 BASRTL.SRCJBASMATMUL .MAR; 1 (5) 
$328 913 FLOAT_TO_FLOAT: SBASSMAT_MUL F, F 
749° 914 


J 12 
ASSMAT_MUL §aBASHC matrix multipt 15-SEP=1 
ae SRAT MO = Mul ttiply 2 arrays giving 6-SEP-19 
fF, 


3243 816 FLOAT_T0_DOUBL: SBASSMAT_MUL 


Be 1:50:88 YGASRTL.sReSBASRATRUL.man:1 88° {3 


K 12 
BASIC matrix multiply 15-SEP-1984 23:4 
BASSHAT : Bt 9330 


~MUL 


= Multiply 2 arrays giving 6-SEP-1984 10:3 


3970 919 FLOAT_T0_GFLOA: SBASSMAT_ MUL F, G 
9A 920 


0:28 


AX/VMS Macro V04-00 
BASRTL.SRC JBASMATMUL .MAR; 1 


L 12 


adn GaSSRAT_RUL''* Rultiply 2 arrays giving '-SEP=1986 10:80:08 UBASRTL.SRESBASMATAUL.man:1 %* {8 
pre 988 FLOAT_TO_HFLOA: SBASSMAT_MUL F, H 


M12 


A T_MUL BAS! matrix multiply i -SEP-1 75 AX/VMS Macro v04-00 Pa 47 
ta baSShay put Mult iD ty 2 arrays giving 6-SEP-19 38 9; $5 £38 EBASRTL. SRE IB JBASMATMUL .MAR; 1 ” (5) 
DC4 5 i+ 
71 : $i. i; Sourcel array is a double array. Now differentiate on the source2 type. 
De, 958 
05 06 02 A3 F 3004 9 § ROUBLE: CASEB paces OTD By Teel WOSCSK_DTYPE_B, #<DSCSK_DTYPE_D = DSCSK_DTYPE_B> 
8 "3009 9 ? 1$: «WORD DOUB _BYTE- Is 3 code for byté dtype 
‘ 3008 9 oW DOUBL “WORD=1 ; code for word dtype 
477° 3000 9 § ‘ DOUBL -10- LONG-1$ ; code for long dtype 
43 DCF 69 oW ERR_DATTYPERR-1$ ¢ quad not supported 
69C' 3DD1 934 " DOUBLE_TO_FLOA-1$ 3 code for float dtype 
08ci' bbs 4 5 «WORD DOUBLE-TO_DOUBL-1$ 3; code for double dtype 
DD 338 3+ 
+44 338 ; G and H floating fall outside the range of the CASEB. 
Dd5 «= 940 * 
18 02 Ay 9} +4 33) cure DSCSB_DTYPE(R3), #DSCSK_DTYPE_G 
OAD8 31 +4 | 348 BRW DOUBLE _TO_GFLOA 
1c 02 as 9 DDE 345 2$ cups DSC$B_DTYPE(R3), #DSC$K_DTYPE_H 
0p02 31 + 4 BRW DOUBLE _TO_HFLOA 
18 02 A 7} 4 rH 3$ Hag pscse. DTYPE(R3), #DSCS$K_DTYPE_DSC 
53 04 A353 DO 3DED 951 MOVL 4(R3), R3 3 R3 <== addr of descriptor 
D1 30F 93¢ BRB DOUB BLE : CASE again on dtype in desc 
C31F 31 3DF3 gee 48 4$: BRW ERR_DATTYPERR 
3re 956 i Now type of sourcel and source2 arrays are known. Use the macro to 
th gee i 3 _generate the code for each case 


N 12 


ASSMAT_MUL BASIC matrix multipl 15-SEP-1984 7:5 AX/VMS Macro v04-00 Pa 4 
35 ~ BASSMA MUL - Muttinty 2 arrays giving 6-SEP-19 3B 9; 30; 38 YBASRTL. SRC JBASMATMUL .MAR; 1 - ‘8, 
3DF6 960 DOUBLE_TO_BYTE: SBASSMAT_MUL dD. B 
4018 961 


ec aa ee Oe Co ee ee ee eee 


B 13 
ASSMAT_MUL BASIC matrix multipl 15-SEP-1984 23: AX/VMS Macro Vv04-00 P 49 
pass BASSMAT MUL = Mult ipty 2 arrays giving 578Eb= 138 ; ‘46: 38 _— | 


4018 963 DOUBLE_TO_WORD: SBASSMAT_MUL D, W 
4240 964 


0: BASRTL. SRE TBASMATMUL. MAR; 1 (5) 


— 
ay 


BASSMAT 
424 966 DOUBLE_TO_LONG: SBASSMAT_MUL 


C 13 
BASIC inatrix multiply 
UL Multiply 2 arrays giving 


446 


15-SEP-19 
6-SEP=-19 


4 
4 


De L 


i 


3: 


0 


AX/VMS Macro v04-00 


40:33 EBASRTL. SRCIBA: MATMUL MAR; 1 


Page 


5 
( 


0 
5 


) 


nem 


TLMUL 


“esa 


D 13 
BASIC matrix putt L 15-SEP 
ba SSMAT 4s 6-SEP 


- Multiply 2 arrays giving 
38° DOUBLE _TO_FLOA: SBASSMAT_MUL 


=198¢ 10:40 


D, F 


28 


AX/VMS Macro V04-00 
BASRTL. SREIGASMATRUL. MAR; 1 


Page 


5 
( 


1 
5 


) 


E 13 
BASIC matrix multipl 15=SEP=1 4 23: 36738 AX/VMS Macro v04-00 Page 52 
BASSMAT _MUL = Mult tiely 2 arrays giving 6=SEP=1984 10:30:23 ([BASRTL.SRCJBASMATMUL MAR; 1 (5) 
reo 444 DOUBLE_TO_DOUBL: SBASSMAT_MUL D. D 
4886 97 


F 13 
| BASSMAT_MUL BASIC matrix multiply 15-SEP= 
1-021 BASSMAT 


1 
MUL = Multiply 2 arrays giving 6-SEP-1 
4886 975 DOUBLE_TO_GFLOA: SBASSMAT_MUL D 
4AEQ =: 9776 


:50 VAX/VMS Macro v04-00 Page 53 
3 (CBASRTL.SRCJBASMATMUL .MAR; 1 (5) 


——_—-~ 


G 13 
BASIC matrix multipl 15-SEP-19 


BASSMAT T_MUL 
4AED 


- Muttinty 2 arrays giving 6-SEP-1984 


Bt 10:80:28 EBASRTL. SRC JBASRATAUL MAR: 


978 DOUBLE_TO_HFLOA: SBASSMAT_MUL D, H 


Page 


54 
(5) 


BASSMAT _MUL 
eat” 


05 06 


18 
1¢ 


18 
53 


02 A3 : 


02 A 91 
0 12 
OAEF 31 
02 a3 91 
03 12 
Op11 = 31 
02 a3 91 
06 is 
04 A3 f' 


@ 
w 
_—J 
o 
w 
— 


z2 


SELOAT: CASEB 
$: i) 


a PPL LLL ELSES Pe 
eer re rnerrre | 
| 


J 
COLL LE EMOMO SP OWN #3 
= 


—— $$ -$DOOCCSCOCCCOoOwovNn 
WO CONOULS WhO 0e 


COOOCOOCOOOoOoOoOoOooowowono 
Se teee 


et es a es 


ti 
tiply 


 Sourcel array is a gfloat array. Now differentiate on the source2 type. 


U 2 arrays giving .* 6-S Epa] 3b 9: $6: 38 YBASRIC. § L.S RE SBASMATRUL MAR: 1 (5) 


H 13 
Page 55 


#OSCSK_DTYPE_8, _M<DSCSK _DTYPE_D = DSCSK_DTYPE_B> 
code for byté dtype 

code for word dtype 

code for long dtype 

quad not supported 

code for float dtype 

code for double dtype 


3% 
a$, ; G and H floating fall outside the range of the CASEB. 


DSCSB_DTYPE(R3), #DSCSK_DTYPE_G 


GFLOAT_TO_GFLOA 

DSC$B_DTYPE(R3), #DSCS$K_DTYPE_H 

GFLOAT_TO_HFLOA 

DSC$B gue #DSCSK_DTYPE_DSC 

4(R3) 3 R3_<-- addr of descriptor 
GFLOAT : CASE again on dtype in desc 


ERR_DATTYPERR 


Now type of sourcel and source2 arrays are known. Use the macro to 
;_generate the code for each case 


I 13 
A T ASI tri Ltiply SEP=1984 23:47:5 AX/VMS Macro v04-00 Pa 5 
eo va ba Beaal wu : Multiply ly 2 arrays giving aes 71 9; :30: £33 EBASRTL. SRCIBASMATMUL MAR; 1 ne 28) 


4D45 1015 GFLOAT_TO_BYTE: SBASSMAT_MUL —G, B 
4F6E 1016 


113 
matrix multipl 15-SEP-1984 7:5 AX/VMS Macro v04-00 57 
MUL’ = Multiply 2 arrays giving ‘G-seP-1984 t TB80:38 POKES TRESS ORT manss OO 2; 


$5 1018 GFLOAT_T0_WORD: SBASSMAT.MUL — G, W 


K 13 
; BASIC matrix multiply 15*SEP=1984 23:47:50 VAX/VMS Macro v04-00 Page 58 
BASSMAT_MUL = Multiply 2 arrays giving B-SEP 19 4 10:30:2 CBASRTL.SRCJBASMATMUL .MAR; 1 (5) 
5197 1021 GFLOAT_TO_LONG: SBASSMAT_MUL G, L 
53C0 1022 


egg 


L 13 
BAS! matrix multipl 15-SEP 
BASSHA MUL Muttipty 2 arrays giving 6-S 


$300 1024 GFLOAT_T0_FLOAT:$BASSMAT MUL 


p-198¢ 10:80:28 EBASATL. SR 


G, F 


v04=-00 
JBASMATMUL MAR; 1 


Pa 59 
” (5) 


wart ba BeAAT RG parr ts multiply shi = 3 P= aa 9 1: $i: 3 y 


= Multiply 2 arrays giving 
7 GFLOAT_TO_DOUBL : $BASSMAT_MUL 
Bic 1055 PtOMT 10. : 


AX/VMS Ma 
BASRTL.SR 


ar 


v04-00 
SMATMUL .MAR; 1 


roe 


1 
$MAT_MUL - Multiply 2 arrays giving 6-SEP-1984 CBASRTL.SRCJBASMATMUL .MAR; 1 5) 


581C 1030 GFLOAT_TO_GFLOA:$BASSMAT_MUL G, G 
SA4? 1031 


N 13 
BASSMAT_MUL BASIC matrix multiply 15-SEP-1984 #3397330 VAX/VMS Macro v04-00 Page 6 
| bas 0:30:2 ( 


egg 


i best 
ASSMAT _MUL 


fa) 


B 14 
C matrix multiply 15-SEP-19 


= Multiply 2 arrays giving 6-SEP-1984 


193? GFLOAT_TO_HFLOA:S$BASSMAT_MUL G, H 


Bs 70:30:25 


AX/VMS 
BASRTL. SRE 


0 V04-00 
JBASMATMUL MAR; 1 


ne 


C 14 


BASSMAT_MUL ; BASIC matrix multipl 15-SEP-1984 23:47:5 AX/VMS Macro V04-00 Page 63 
3-021 - BASSRAT_MUL’ = MULTIDLY 2 arrays giving ‘e-seP-19e¢ 10:30:29 FBASRTL SAC IBASHATAOL maR:1 2% $3 
c78 1 3+ 
l i ; Sourcel array is an hfloat array. Now differentiate on the source2 type. 
C78 1039 ° 
05 06 02 A3 F 5C7 1940 HFLOAT: CASEB DSCSB_DTYPE(R3), MDSCSK_DTYPE_B, #<DSCSK_DTYPE_D = DSCSK_DTYPE_B> 
09 d' 5¢7D 104 1$: <WORD HFLOAT_TO_BYTE-{$ ; code for byte dtype 
6° SC7F 1 ts «WORD HFLOAT_TO_WORD-1$ ; code for word dtype 
47F* 5(81 13 ° HFLOAT_ TO_LONG-1$ ; code for long dtype 
AGOR 583 1044 «WORD ERR_DATTYPERR-1$ : quad not supported 
06A8' cee oF «WORD HFLOAT_TO_FLOAT-1$ ; code for float dtype 
08D1* 5¢87 1 $ «WORD HFLOAT_TO_DOUBL-1$ ; code for double dtype 
at 89 1948 -+ 
S44 1049 : G and H floating fall outside the range of the CASEB. 
gS 1838 © 
1B 02 a3 } 589 1936 we DSCSB_DTYPE(R3), #DSCSK_DTYPE_G 
1 3C8F 1054 BRW HFLOAT_TO_GFLOA 
A } s¢90 1928 2s coe DSCSB_DTYPE(R3), MDSCSK_DTYPE_H 
op0B «31 398 1058 BRW HFLOAT_TO_HFLOA 
18 02 aa } ocoe 1399 3$ ee PECMS PT WE ED, #DSCSK_DTYPE_DSC 
53 04 AS Q 5CAl 1996 MOVL 4(R3), R3 ; R3 <== addr of descriptor 
D1 1 ace ier BRB HFLOAT ; CASE again on dtype in desc 
A46BOU331 ap 1302 4$: BRW ERR_DATTYPERR 
5SCAA 1067 : Now type of sourcel and source2 arrays are known. Use the macro to 
pe 1368 3; generate the code for each case 


ah | 
o ' 
oO 
nm Oo 
od 
> m 
ww 
-c-0 WO WO 


D 14 
ASSMAT_MUL ; BASIC matrix multipl 15-SEP-1984 23:47: AX/VMS Macro V04-00 Pa 64 
passe BASSAAT_RUL’ = MULEGDLY 2 arrays giving ‘e-sEP-19b6 {0:40.29 FBASNTE eRe TBASMATHoL.man;1 2% $8, 
SCAA 
5ED3 


1071 HFLOAT_TO_BYTE: SBASSMAT_MUL H, B 


egg 


pasts matrix multiply 


BASSMAT 


core 


~MUL 


5:40:23 


i os ad bs 9 


E 
- Mult tiply 2 arrays giving 6-SEP-1984 


O74 HFLOAT.TO_WORD: SBASSMAT MUL 


AX/VMS Macro v04-00 
BASRTL.SRCJBASMATMUL .MAR; 1 


Page 65 


(5). 


———— - — 


— 


F 14 
BASSMAT_MUL BASIC matrix multipl 1 P=19 
1-021 


5-SE : AX/VMS Macro v04-00 P 66 
BASSMAT T_MUL * Rultipty 2 arrays giving 6-SEP-19 : g ree 


BASRTL.SRCJBASMATMUL .MAR; 1 (5). 


@o 
- FS 
=~ 
ow 
we 
on 
mow 


60FC 1077 HFLOAT_TO_LONG: SBASSMAT_MUL +H, 
6325 1078 


—_—_—— 


G 14 
ASSMAT_MUL ; BASIC matrix multiply 15-SEP-1984 23:47:5 AX/VMS Macro v04-00 Pa 67 
at BASSMAT_MUL = Multip ty 2 arrays giving g-$Ep= 198 9; :30: £33 (BASRIL. SR SRC JBASMATRUL .MAR; 1 ~ (5). 
$323 1980 HFLOAT_TO_FLOAT:$BASSMAT_MUL H, F 
654E 1081 


eganernn 


H 14 


GASSAAT MUL” '= RULEIDLY 2 arrays giving '6-SEP-198¢ 10:30:28 EBASRIL.SRCIBASMATMUL.MaR;1 8° (8). 
M.D 


6546 1085 HF LOAT_TO_DOUBL : $BASSMAT_MUL 


1 14 
RA SMAT_MUL BASIC matrix multipl 4° mi P-19 
21 BASSMAT MUL * uttipty 2 arrays giving P=19 
H 


6777 1986 HFLOAT_TO_GFLOA:S$BASSMAT_MUL 
69A6 1087 


49 $3: 47: +30 oben Macro v04-00 Page 69 
4 10:30:2 BASRTL.SRCIJBASMATMUL . MAR; 1 (5). 


AaBASHE patrix multiply ip pee aes 95792 '38 


= Multiply 2 arrays giving 
69A6 1089 HFLOAT_TO_HFLOA:SBASSMAT_MUL 


AX/VMS Mac 
BASRTL.SRC 


ro 
JBA 


v04- 
SMAT 


00 


MUL .MAR; 1 


re 1 


K 14 
BA SMAT_MUL ; BASIC matrix multipl 15-SEP-1984 23: AX/VMS Macro v04-00 we 71 
21 GASSAAT MUL’ = RULEIDLY 2 arrays giving ‘O-sePol98e TOr20;5% FOASHTS cae ToatonT on mar:1 29° Zh, 
6BD1 1091 i* 
peed ! {ee ; Add has been in byte. Determine destination type to convert to dest. 
6BD1 1094 °- 
68D1 1095 DEST_CASE B: 
50 56'° H 980} ! My we current. sum(R5), RO i get " te —. 9 As 
05 06 02 A6 BF 6BD8 1 a8 5$ CASEB DSC§B_DTYPE(R6), #DSCSK_DTYPE_B, * pebstsx Mis hs aes _D - eee _DTYPE_B> 
1A?" 6BDD 1099 1$ -WORD STORE BYTE-1$ 3 no conversion needed 
ee. 6BDF 1100 «WORD DEST_B_TO_wW-1$ ; code for word dtype 
6BE1 1101 «WORD DEST-B-TO-L-1$ : code for long dtype 
5 pert 11 ¢ «WORD ERR_BATTYPERR-1$ 3 quad not supported 
Be * 6BE5 11 ~WOR DEST_B_TO_F-1$ ; code for float dtype 
580° 6BE7 1104 «WORD DEST_B_TO_D-1$ ; code for double dtype 
8s 
6BE9 11 5 : G6 and H floating fall outside the range of the CASEB. 
ss 188 
18 «(02 %} yo 1119 oe DSCSB_DTYPE(R6), #DSCSK_DTYPE_G 
O6CA 531 eect BL BRwW DEST_B_T0_G 
1¢ 02 ag 3 ear? 1iN3 23: ee DSCSB_DTYPE(R6), #DSCSK_DTYPE_H 
O7DA 31 ear 1118 BRwW DEST_B_TO_H 
18 02 ry4 bf core 1333 3$ eee DSCSB_DTYPE(R6), #DSCSK_DTYPE_DSC 
56 04 A6 6 6C01 1120 MOVL 4(R6), R6 R6 <-- addr of descriptor 
D1 11 oc8e 1% 1 BRB 5$ CASE again for dtype in desc 
9508 31 6C0O7 11 § 4$ BRW ERR_DATTYPERR 
eh 1124 
sok N82 
gene 1 3 : Add has been in word. Determine destination type to convert to dest. 
6COA 1199 ° 
6COA 1130 DEST _CASE vs 
50 os" as 00 ocRA 1 1 nove currens sun(RS), RO : get. ’ otf store ty ne ne 
ina 
05 06 02 A6 F tat 11 : 5$: CA SEB DSC$B_DTYPE(R6), #DSCSK _DTYPE_B, * pedsesk _DTYPE meant: DSCSK _DTYPE_8> 
150° ecg 1134 1$: -WORD DEST@ TO B- ; code for byté dtype 
; 1° 6€18 1135 -WORD STORE BORD™ 18 3 MO conversion needed 
57° 6C1A 11 § « DEST @ TO L-i$ ; code for long dtype 
4FF 6C1C 113 . ERR 7A TIVBERRS1S 3 quad not supported 
Besa. 6C1 1138 -WORD DEST_W_TG_F-1$ : code for float dtype 
550° 6C 1139 WORD DEST“W>TO- D-1$ 3; code for double dtvee 
$53 1141: 
ot 1106 3 > G and H floating fall outside the range of the CASEB. 
6C 1144 * 
18 02 Ag ? of 136? ey RSCSB_DTYPE (RO). #DSCSK_DTYPE_G 
0697 4 6C 1143 BRw DEST_W_TO_G 


a 


L 14 
ASSMAT_MUL BASIC matrix multipl 7 SEP- AX/VMS o v04-00 e 7 
35h BASSAAT MUL’ = MULEIOLY 2 arrays giving ‘O-sEPa19R6 {0:90.29 FOASNTr eae TBASMATHoL.man;1 2% 7, 
C2B 114 
1c 02 03 % e¢ ? 1148 2$: me SCSB_DTYPE(R6), #DSCSK_DTYPE_H 
O7A if of } Hi BRw DEST_W_TO_H 
| 18 02 ap 93 6¢ , 1188 3$ cups DSC$B_DTYPE(R6), #DSCSK_DTYPE_DSC 
56 04 A6 8 6C3A 1155 MOVL 4(R6), R6 3; R6 <-- addr of descriptor 
vi 11 oz 1128 BRB 5$ ; CASE again for dtype in desc 
94D2 31 + 6C4 1128 4$: BRW ERR_DATTYPERR 
gte8 1100 s+ 
ote 1188 3 Add has been in long. Determine destination type to convert to dest. 
6643 1163 ° 
6C43 1164 DEST_CASE_L: 
50 os" oH 4 otés 1393 nov. gycrens sum(R5) RO : get. a to's were, f9 Az 
e or 
05 06 02 A6’ 8F 6C4A 118? 5$: CASEB osc $s DTYPE(R6), #DSCSK_DTYPE_B, “#<DSC8K cof ree oP °° ate _DTYPE_B> 
O11C* 6C4F 1168 18: WORD DEST_C_T0_B-1$ 3 code for byte dtype 
pete. 6C51 1169 WOR DEST L_ TO We ; code for word dtype 
0338' 6C€53 1170 WORD STORE_CONG-1 3 Mo conversion needed 
94€6 6C55 1171 «WORD ERR_DATTYPERR-1$ 3 quad not supported 
eth 657 11%¢ WORD DEST_L_TO_F-1 3; code for float dtype 
052D° 4TH Ma - -WORD DEST_L-T0_D-1$ 3; code for double dive 
6C5B 1175 ; G and H floating fall outside the range of the CASEB. 
step ine 
18 02 A6 91 6C5B 1128 CMPB DSCSB_DTYPE(R6), #DSCSK_DTYPE_G 
03 12 6CSF 117 BNEQ 
0664 31 ore} 1189 BRW DEST_L_TO_G 
1¢ 02 AB «91-6064 118¢ 2$: CMPB SCSB_DTYPE(R6), #DSCSK_DTYPE_H 
03 12 6C68 118 BNEQ 
0774 31 eo ry 1138 BRW DEST_L_TO_H 
18 02 ae 7 ocee 1138 3$: cee DSCSB_DTYPE(R6), #DSCSK_DTYPE_DSC 
56 04 A6 46 6C73 1188 MOVL 4(R6), R6 ; R6 <== addr of descriptor 
ee | ocr? 1133 BRB 5$ ; CASE again for dtype in desc 
9499 31 6C79 1191 4$: BRW ERR_DATTYPERR 
6C7C 1336 
sre Ne 
ecre 1186 : Add has been in float. Determine destination type to convert to dest. 
6c7C 1197 ° 
6C7¢ 1198 DEST_CASE_F: 
50 1¢ aS 50 6C7C 119 MOVE current_sum(R5), RO : get # to yore. ty AF 
56 =SA ” 6C 1 $9 MOVL R10, R6™ ae ies 
05 06 02 A6 F 6C 1201 5$: CASEB osc $B DTYPE(R6), #DSCSK_DTYPE_B, “#<DSCSK DIYPE_D =- eS _DTYPE _B> 
oes 6C 1 8 1$: WORD DEST_F_TO_B-1$ : code for byté dtype 
1E7° +3 A 1 WORD DEST"F"-TO-W-1$ : code for word dtype 
2EA* 6C8C 1204 «WORD DEST_F_TO_L-1$ ; code for long dtype 
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BASSMAT_MUL ; BASIC matrix multipl 15-SEP-1984 23:47:5 AX/VMS Macro v04-00 P 
wits} BASSMAT_MUL - Multiply 2 arrays giving 6-SEP-1984 93:40:38 EBASRTL. SREIGASMATMUL MAR: 1 a 
948D 6C 1205 «WORD ERR_DATTYPERR-1$ 3 quad not supported 
Be 1° 3 1 § «WORD STORE_FLOAT-1$ : bm senvers ton needed 
FD 34 : . «WORD DEST_F_TO_D-1$ ; code for double dtype 
| oboe ! % 3 G and H floating fall outside the range of the CASEB. 
6(94 1211 ° 
18 02 ag bf oboe ! \¢ oars * a ahthaaaaat #DSCSK_DTYPE_G 
0681 if ora ! 13 BRW DEST_F_TO_G 
1¢ 02 " + tH : 1g 2$: coe QgCSO_OTVPE CAE). #DSCSK_DTYPE_H 
0741 31 ork? ! 1% BRW DEST_F_TO_H 
18 02 hve 33 TeAe : Y 3$ care DSCSB_DTYPE(R6), #DSCSK_DTYPE_DSC 
56 04 A6 00 6CAC 1 ¢ MOVL 4(R6), R6 3; R6 <-- addr of descriptor 
D1 11 th : ? BRB 5$ ; CASE again for dtype in desc 
9460 31 ete 1225 4$ BRW ERR_DATTYPERR 
6CB5 1226 
Seas 1598 + 
ote : 4 ; Add has been in double. Determine destination type to convert to dest. 
6CB5 1231 ° 
6cBs 1 ; DEST_CASE_D: 
50 1C AS 70 6CBS 123 mOvVD current_sum(R5), RO i; get # to store in RO&R1 
SA OC AC DO 6CB9 1234 MOVL dest_matrix(AP), R10 3 point to dest matrix 
56 SA 00 6CBD \$32 MOVL R10 3 Save original pointer 
05 06 02 A6 BF 6CCO 1236 5S: CASEB DSC$B_DTYPE(R6), #DSCSK_DTYPE_B, #<DSCSK_DTYPE_D - DSC$K_DTYPE_B> 
0080" 6CC5 \s37 18: «WORD DEST_B_TO_B-1$ : code for byte dtype 
OIAF* 6CC? 1238 «WORD DEST_D_TO_W-1$ :; code for word dtype 
oc 8s 6cc9 1239 0 DEST D TO L-1$ ; code for long dtype 
9450 6CCB 1240 -WORD ERR_BATTYPERR-1$ 3 quad not supported 
0385' 6CCD 1241 «WORD DEST_D_TO F-1 ; code for float dtype 
0516° occ : 4s Be -WORD STORE_BOUBLE-1$ : no conversion needed 
tt \ rf ; G and H floating fall outside the range of the CASEB. 
6CD1 1246 * 
1B 2 a8 7 eco! { ¢} Hae DSCSB_DTYPE(R6), #DSCSK_DTYPE_G 
OSFA if eit } i BRW DEST_D_TO_G 
1C 02 ag | 7508 : 2) 2s ce QSCEO_DTVPE (RO), #DSCSK_DTYPE_H 
070A 4 a : 38 BRWwW DEST_D_TO_H 
18 02 a6 9} ece$ \ 22 3$ ee DSCSB_DTYPE(R6), #DSCSK_DTYPE_DSC 
56 04 A6 DO 6CED 1 3§ MOVL 4(R6), R6 3 R6 <-- addr of descriptor 
D1 11 oceD ' 38 BRB 5$ ; CASE again for dtype in desc 
9423 31 6CEF 1 99 4$ BRW ERR_DATTYPERR 
6CF2 126 
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mu mutt ip ty 2 arrays giving % Fa Pa] 3b 19; 46: 39 YOASRTL se 


o V04-00 Page 74 
SRE SBASMATMUL. MAR; 1 (6) 


1 6¢ 
1265 5+ 
! be ; Add has been in gfloat. Determine destination type to convert to dest. 
1266 * 
1 $9 DEST_CASE_G: 
1268 MOVG current_sum(R5), RO : get - te store 19 Ad 
1 4 OVL R10, R6 ave rAginal poi 
1270 5$: CASEB DSC $B DTYPE(R6), #DSCSK_DTYPE_B, "#<OSC8K D ry o- DSe$K _DTYPE_8> 
1271 1$: «WORD DEST_G_TO_B-1$ : code for byté dtype 
1 es -WORD DEST_G_TO_W-1$ 3 code for word dtype 
127 «WORD DEST GTO L-1$ ; code for long dtype 
1274 -WORD ERR_BATTYPERR-1$ 3 quad not supported 
1275 «WORD DEST_G_TO_F-1 i code for float dtype 
: 6 «WORD DEST_G_TO_D-1 :; code for double dtype 
1278 ;+ 
: ; G and H floating fall outside the range of the CASEB. 
1281 © | 
1 82 CMPB =: DSC$B_DTYPE(R6), #DSCSK_DTYPE_G | 
128 BNEQ 2$ 
84 BRW STORE_GFLOAT | 
1386 28 CMPB -«DSCSB_DTYPE(R6), #DSCSK_DTYPE_H | 
1287 BNEQ 3$ 
\so8 BRW DEST_G_TO_H 
1290 3$ CMPB = DSC$B_DTYPE(R6), #DSCSK_DTYPE_DSC 
1291 BNEQ 4$ | 
1 38 MOVL 4(R6), R6 3 R6 <-- addr of descriptor 
Z again for ype in desc 
: + BRB 5$ CASE for dt d 
1095 4$: BRW ERR_DATTYPERR 
1599 | 
3+ 
3 as been in Oat. Determine destination type to convert to dest. | 
' 4 Add has b in hfloat. Determine destination t t t to dest 
1300 ° | 
1301 DEST_CASE_H: 
1 8 MOVH current_sum(R5), RO : get . te store in yd 
130 MOVL R10, R6 f hgine point 
1304 5$ CASEB SC $B DTYPE(R6), #DSCSK_DTYPE_B, *weDseSk _DTYPE ‘ - DSe$k _DTYPE_B> 
1305 1$ WORD DEST_A_TO_B-1$ ; code for byté dtype 
1 $ «WORD DEST_H_TO_W-1$ ; code for word dtype 
1 -WORD DEST H_ TO L-1 ; code for long dtype 
1308 -WORD ERR_SATTYPERR-1$ 3 quad not supported 
1309 WORD DEST_H_TO_F-1 ; code for float dtype 
: + -WO DEST_H_T0_D-1 ; code for double dtype 
1312 ;+ 
3: G an coating fall outside the range o 2 . 
: 18 G and H floati fall outside th f the CASEB 
1315 * 
1 1 CMPB DSCSB_DTYPE(R6), #DSCSK_DTYPE_G 
131 BNEQ 
1318 BRW DEST_H_T10_G 


BASSMAT_MUL 
i981 = 
1¢ 
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; BASIC matrix multipl 
BASSMAT_MUL = Multiply 
6D4E 1319 
02 a6 91 gD4e 1320 2$ CMPB 
0 \¢ 6052 1321 BNEQ 
ape pak ES Sy Aan 
02 A6 91 «+6057 «1324 3S CMPB 
06 if 6058 1325 BNEQ 
04 A6 DO 605d 1 : MOVL 
Di «11 «6061 «1 BRB 
6065 1328 
Q3AF 31 6063 1329 4S: BRW 
Spee 1331 DEST_w.10.B 
50 33 6066 1 & -* COTwB 
19 11 6069 1 BRB 
388 1H vesrt.to. 
50 F6 6D6B 1 36 -"" COTLB 
14 11 6D6E 133 BRB 
$078 1339 DEST_F_TO_B: 
50 48 6070 1340 corr FB 
OF 11 6073 1341 
6D75 1348 
6075 134 ie B: 
50 68 6075 1344 CTD DB 
OA 11 6078 1345 
6D7A 1346 
6D7A 1347 DEST_G_TO_B: 
50 48FD 6D7A 1348 C0TGB 
04 11 6D7E 1349 BRB 
6080 1350 
6D80 1351 DEST_H_TO_B: 
50 68FD 6D80 1358 COTHB 
6084 135 
6084 1354 
6D84 1355 STORE_BYTE: 
SA DO 6D84 1356 MOVL 
18 AS DO 6D87 1357 MOVL 
1405 DO 6D8B 1358 MOVL 
50 90 6D8F 1359 MOVB 
6093 1360 STORE 
05 6664 1361 RSB 
6£65 136¢ 
6£65 1363 DEST_B_TO_w: 
50 99 6665 1364 COTBwW 
1D 11 6668 1365 BRB 
6E6A 1366 
6E6A 1367 DEST_L_TO_w: 
50 F7 6E6A 1368 COTLW 
18 11 66D 1369 BRB 
6E6F 1370 
6E6F 1371 DEST_F_TO_w: 
50 49 6E6F 1 es COTFW 
1 11 6€72 137 BRB 
6674 1374 
6£74 1375 DEST_D_TO_w: 


2 arrays giving 
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sc #DSCSK_DTYPE_H 
TORE _HFLOAT 
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5 
D 
4 
G(R 
5$ 
ERR_DATTYPERR 
RO, RO 
STORE_BYTE 
RO, RO 
STORE_BYTE 
RO, RO 
STORE _BYTE 
RO, RO 
STORE _BYTE 


RO, RO 


RO, RO 
STORE _WORD 
RO, RO 
STORE_WORD 


RO, RO 
STORE _WORD 
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SBASMATMUL.MAR:1 °° 


R6 <== addr of descriptor 
CASE again for dtype in desc 


convert 
go store 


2 convert 
» go store 


+ convert 
: go store 


i convert 
3 go store 


convert 
qo store 


conve 
fall inte store 


pointer to dest descriptor 
current row 
current column 


go continue loop 


convert 
go store 


convert 
go store 


convert 
go store 
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BASSMAT_MUL ; BASIC matrix multipl 15-SEP-1984 23:47:5 AX/VMS Macro v04-00 Pa 76 { 
4 21 BASSMAT_MUL - Multiply 2 arrays giving 3-8Ep-y ope 98: 40.3% EBASRTL. SREIBASMATHOL..MAR: 1 ” (6) | 
50 2C AS 66 6E74 137 DIvD2 scale(R5), RO : descale for dest 
+ 69 6€78 1 5 CVTDW RO, R 3 convert to word 
A 11 6€7B 1 £8 BRB STORE_WORD ¢ go store 
6€70 137 
6E7D 1380 DEST_G_TO_w: 
50 50 49FD 6E7D 1381 CVTGw = RO, RO 3 convert 
6 «11 at 86 BRB STORE _WORD ; go store 
6E83 1384 DEST_H_TO_w: 
50 50 69FD 6€ 1385 CVTHW RO, RO 3 convert 
6E 1386 3; fall into store 
6€87 1 iH STORE _WORD: 
51 SA DO 6687 1388 MOVL R10, R1 3; pointer to dest descriptor 
2g 18 AS DO 6E8A 1389 MOVL current_i(R5), Re 3 current row 
5 14 aS DO oEee 1390 MOVL current 4 SE? R 3 current column 
GA AE 50 BO 6€9 1391 MOVW RO, DATA(SP) 
ore 1 5 STORE WwW ; store — 
05 6F67 139 RSB 3 go continue ioop 
6F68 1394 
6F68 1395 DEST_B_TO_L: 
50 50 98 6F68 1396 CVTBL = RO, RO 3; convert 
1D 11 6F6B 1397 BRB STORE _LONG : go store 
6F 6D 1338 
6F6D 1399 DEST_W_TO_L: 
50 50 32 6F6D 1400 CVTwL 3=RO, RO : convert 
18 11 6F70 1401 BRB STORE _LONG 3 go store 
oF es iene 
6F72 1403 DEST_F_TO_L: 
50 50 4A 6F72 1404 CVTFL RO, RO 3 convert 
15. 11 «6F75) «(1405 BRB STORE_LONG 3; go store 
6F77 1406 
6F77 1407 DEST_D_TO_L: 
50 2C AS 66 6F77 1408 DIVD2 scale (R5), RO ; descale for dest 
0 50 6A 6F7B 1409 CVTOL §«6—RO, R 3: convert 
OA 11 6F7E 1410 BRB STORE _LONG ; go store 
6F80 1411 
6F 80 Vals DEST_G_TO_L: 
50 50 4AFD 6F80 141 CVTGL RO, RO 3 convert 
04 11 GFB4 1414 BRB STORE LONG ; go store 
6F86 1415 
6F86 1416 DEST_H_TO_L: 
50 50 6AFD 6F86 1417 CVTHL RO, RO 3 convert 
6F BA 1213 ; fall into store 
A oy 141 
6FBA 1420 STORE_LONG: 
51 SA DO 6F8A 1421 MOVL R10, R1 3; pointer to dest descriptor 
28 18 AS) DO 6F8D 14 ¢ MOVL current_i(R5), Re 3 current row 
5 14 AS DO 6F91 14 MOVL current JiR). R 3 current column 
4A AE 50 DO 6F95 1424 MOVL RO, DATACSP) 
44 1425 STORE L ; store | 
05 706A 14 § RSB 3 go continue Loop 
706B 14 
7068 1428 DEST_B_TO_F: 
50 50 4C 7068 1429 CVTBF «RO, RO 3: convert 
7. i PF 6 1430 BRB STORE_FLOAT : go store 
70 1431 
7070 14352 DEST_W_TO_F: 
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BASSMAT_MUL ; BASIC matrix multipl 15-SEP-1984 23:47:5 AX/VMS Macro v04-00 Pa 77 
wth BASSMAT_MUL = Multiply 2 arrays giving 6-SEP- 71984 $9:40:39 BASRTL.SRCIBASMATMUL.MAR;1 > (6) | 
50 50 4D 7070 1433 CVTWF 8 6RO, RO 3 convert 
4 «11 f : 14 4 BRB STORE_FLOAT ; go store 
7075 1436 DEST_L_TO-F: 
50 50 4€ 7075 14 $ “" COTLE RO, RO 3 convert 
OF 11 7078 1438 BRB STORE_FLOAT ; go store 
, th 1220 DEST_D_TO_F 
36 DM aN 1441 “" cVTOF RO, RO 3 convert 
OA 11° 7070 1446 RB STORE_FLOAT : go store 
707F 144 
707F 1444 DEST_G_TO_F: 
50 50 33FD 7O7F 1445 COTGF RO, RO 3; convert 
04 11 «+7083 1446 BRB STORE_FLOAT 3 go store 
7085 144 
7085 1448 DEST_H_TO _F: 
50 50 F6FD 7085 1449 CUTHF RO, RO 3: convert 
7089 1450 ; fall into store 
7089 1451 
7089 1926 STORE _FLOAT: 
51 5A DO 7089 145 MOVL R10, R1 3 pointer to dest descriptor 
2g 18 AS DO 708C 1454 MOVL current_i(R5), ~ 3; current row 
5 14 AS 00 7090 1455 MOVL cerrent ht R 3 current column 
4A AE 50 50 7094 1456 MOVF 0, DATACSP) 
7098 1457 STORE f 3; store 
05 7169 1458 RSB 3 go continue Loop 
716A 1459 
716A 1460 DEST_B_TO_D 
50 50 6C 716A 1461 C0TB BD RO, RO : save double 
50 2C AS 64 716D 1006 MULD2 = scale(R5), RO : scale for dest - 
7171 +146 3 no integerize necessary 
68 11 7171 1464 BRB STORE _DOUBLE 3 go store 
7173 «1465 
ate: 1466 DEST_W_TO_D: 
0 6D 7173 1467 C0TwD RO, RO 3; save double 
50 2C AS 64 7176 1468 MULD2 = scale(R5), RO ; scale for dest 
717A 1469 3; no integerize necessary 
SF 11 717A 1470 BRB STORE DOUBLE ; go store 
717C «(1471 
717C¢ 1076 DEST_L_TO_D: 
5 50 6€ 717C 147 CUTLD RO, RO 3; save double 
50 2C AS 64 7I7F 1474 MULD2 scale(R5), RO ; scale for dest 
aH 1475 3 no integerize necessary 
56 11 me 1676 BRB STORE _DOUBLE ; go store 
7185 147 
£185 1478 DEST_F_TO_D: 
50 50 56 7185 1479 CUTFD = =6RO, RO : save double 
50 c AS ¢ 7188 1480 MULD2  scale(R5), Rg ; scale for dest 
08 C AS 1 718C 1481 CMPD scale(R5), : is the scale 0? 
® ww Fy 14 $ BEQL 1 ; yes, do not integerize 
00000000'GF 16 7192 14 JSB G“MTHSDINT_R4 3 no, integerize 
41 11 7198 1484 1$: BRB STORE _DOUBCE 3 go store 
719A 1485 
719A 14 $ DEST_G_TO_D: 
719A 14 37 ‘ ; 
719A 1488 3; Note the intermediate conversion to hfloat. 
719A 1489 3° 
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BASSMAT_MUL ; BASIC matrix multipl 15-SEP-1984 23:47: AX/VMS Macro V04-00 Page 78 
ety BASSMAT_MUL - Multiply 2 arrays giving Pat eet 9 3:45:35 YBASRTL. SRE IBASMATHUL.MAR: 1 . (6) 
7 D 19A 1490 MOVL R2, -(SP) 3 save regs which CVTGH 
re D r D 1491 MOVL R35, =(SP) : will eestrey 
S6FD 71A0 1436 CVTGH RO, Re 3 cvt gfloat to hfloat 
QO F7FD 71AQ 149 CVTHD R 3 cvt to desired double 
E DO 71A8 149% MOVL  (SP)+, R3 ; restore regs 
3 DO 71AB 1495 MOVL (SP)+, R2 
50 3¢ A $4 c1AE 1496 MULD2 scale(r5), RO : scale for dest 
08 C Ag 1 71B2 149 CMPD scale(R5), #1 3; scale 0? 
¢ 13° 7186 1498 BEQL STORE BLE ; yes, don't integerize 
7E 4 00 71B8 1499 MOVL R4, -T 3 save R4 
00000000'GF 16 71BB 1500 JSB G*MTHSDINT_R4 : integerize 
4  8€ 0 71C1 138) MOVL (SP)+, R 3: restore R4 
0014 1 #716 15 § BRw STORE _DOUBLE 
71C7.—- 150 
71C7 1504 DEST_H_TO_D: 
50 SO F7FD 71C7 1505 CVTHD RO, RO 3 save double 
50 CAS 64 71CB 1208 MULD2 = scale(R5), RO ; scale for dest 
08 CAS 71 7ICE 150 CMPD scale(R5), #1 ; is the scale 0? 
06 13 7103 1508 BEQL STORE DOUBLE ; yes, do not integerize 
00000000'GF 16 71D5 1509 JSB G*MTHSDINT_R4 3; no, integerize 
7108 1510 3 fall into store 
71DB 1511 STORE_DOUBLE: 
52 SA DO 71DB 1316 MOVL R10, R2 3; pointer to dest descriptor 
53 18 AS 00 pips 151 MOVL current_i(R5), R3 3 current row 
54 14 A5 D0 71E2 1514 MOVL current Bh ae R4 3; current column 
GA AE 50 70 71E6 1515 MOVD RO, DATA(SP) 
71EA 1318 STORE D 3; store 
05 7288 151 $B 3 go continue loop 
72BC 1518 
72BC 1519 DEST_B_T0 G: 
50 50 4CFD 72BC 1360 CVTBG) = RO, RO 3 convert 
2D 11 rset 13s! BRB STORE _GFLOAT ; go store 
73t3 15 : DEST_W_TO_G: 
50 50 4DFD 72C2 1524 CVTwG = =RO, RO 3 convert 
e7 11 72C6 1525 BRB STORE _GFLOAT 3 go store 
72C8 1526 
rec 1527 DEST_L_T0 G: 
50 50 4EFD 72C8 1528 CVTLG RO, RO 3 convert 
21 11. 72CC 136? BRB STORE _GFLOAT ; go store 
rece 1530 
72CE 1531 DEST_F_TO.G: 
50 50 99FD 72C 1338 CVTFG RO, RO : convert 
18 11 72D 1 BRB STORE _GFLOAT ; go store 
7204 1534 
7204 1535 DEST_D_TO_G: 
7204 1536 ;% : ; 
f + 3 i 3 Note the intermediate conversion to hfloat. 
7E 26 DO 7204 1539 MOVL Ree -(SP) 3; save regs which CVTDH 
7E 5 D0 72D7 1540 MOVL R35, -(SP) : 6Uwill greteoy 
50 = 50 sere 72DA 1541 CVTDH = RO, RO ; cvt dbl to hfloat 
29 50 76FD 72D 1366 CVTHG RO, RO 3 cvt to desired gfloat 
5 ge DO 72E2 154 MOVL (SP)+, R3 3 restore regs 
52 E Y 72E 1363 MOVL (SP)+, R2 
0004 1 7268 1545 BRW STORE _GFLOAT 
72EB 1546 
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ASSMAT ; BASIC matri Ltiply os 4 7:5 AX/VMS o V04-00 
beat See: BASSMAT_MUL : Multip ip ty 2 arrays giving 6S Pa} 984 99; 46: 38 LBASRIL. SRC TBASMATMUL. MAR; 1? 
72EB 1547 DEST_H_T0_G 
50 50 76FD 72EB 1 $3 COTHG RO, RO 3 convert 
f +} } rH 3; fall into store 
72EF 1 i STORE_GFLOAT: 
52 SA DO 72EF 1 ¢ MOVL R10, R2 ; pointer to dest descriptor 
7: a2 09 7? ee 1 2 MOVL current_i(R5), R3 3; current row 
4 6A DO 72F6 1554 MOVL current {soe R4 3 current column 
4A AE 50 SOFD f . 132? ores RO, DATA(SP) 
05 7 oS 1387 RSB 3 go continue loop 
5 Be 1238 DEST_B_TO_H: 
50 50 6CFD 7305 1560 svien RO, RO 3 convert 
| ie 1309 1361 STORE _HFLOAT 3 go store 
7308 1308 DEST_W_TO_H: 
50 50 6DFD 730B 1564 CVTWH RO, RO 3 convert 
16 «11 capt 1392 BRB STORE _HFLOAT 3; go store 
73E1 1369 DEST_L_TO_H: 
50 50 6EFD 73E1 1568 COTLH RO, RO 3 convert 
10 #11 Pees 129) BRB STORE _HFLOAT 3 go store 
eae? 1571 DEST_F_TO_H: 
50 50 98FD 73€7 1376 cOTFH RO, RO 3: convert 
OA 11 geen 1377 STORE _HFLOAT ; go store 
73ED 1575 DEST_D_TO_H: 
50 50 32FD 73ED 1576 CVTDH RO, RO 3 convert 
04 11 73F1 1577 BRB STORE _HFLOAT 3 go store 
pare 1358 DEST_G TO_H: 
50 50 56FD oe 1580 =" COTGH RO, RO : convert 
rare 133 : fall into store 
3? 1388 STORE _HFLOAT: 
54 SA DO 73F7 1584 MOVL R10, R4 ; pointer to dest descriptor 
56 14 A5 00 Para 1585 MOVL current_j(R5), R6 : current column 
55 18 AS DO 73FE 1586 MOVL current _1(R5), R5 3: current row 
4A AE 50 70FD oe 138? iene RO, R(SP) 
05 740C 1589 SB 3; go continue loop 
74DD 1590 -END 
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BASSMAT _MUL ; BASIC matrix multipl 15-SEP-1984 23:47:5 AX/VMS Macro V04-00 Pa 83 
Psect synopsis te b-8Ep= 1984 19:45 3 EBASRIL. SR SRC JBASMATMUL .MAR; 1 ot (6) 
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! ;_Psect synopsis : 


Steet te te 


PSECT name Allocation PSECT No. Attributes 
» O88 « 0000000 < 0.) oF 9. NOPIC USR ABS NOSHR NOEXE NORD NOWRT NOVEC BYTE 
SABSS$ 0000000 ( N 
_BASSCODE 00074DD ( 
$ee scene noma aaroccw eso ne oe 
! Performance indicators ! 
Phase Page fau’*s CPU Time Elapsed Time 
Initialization 28 98:00:00 .08 00:00:00.37 
pounene processing 105 0:00:00.61 ss Be BY 
1223 00:00:49.75 00:01:41.26 
Synbol table sort 00:00:02.28 BF 8008-2 
s2 87 00:00:12.69 82:82: 3.90 
eee table output 4 888036 26 0:00:00.59 
Psect synopsis output 4 0:00:00.04 00:00:00.09 
Cross-reference output 0 00:00:00.00 00:00:00.00 
Assembler run totals 2278 00:01:05.71 00:02:23.59 


The vorking set $i tate was 900 pages. 

366875 bytes ( pages) of virtual memory were used to buffer the intermediate code 

There were 70 cae of symbol table peace at locates to hold 479 non-local and 955 Local symbols. 
source lines were read in Pass 1, produci " 91 object records in Pass 2. 

46 pages of virtual memory were used to define 11 macros. 


eect en eee eee eee este sw ee aw ea > 


Macro Library name Macros defined 

-§ SSSDUA . BASRTL.OBJJBASRTL.MLB; 1 g 
$255$DUA28:CSYSLIBISTARLET. MLB; 2 

TOTALS (all libraries) 7 


493 GETS were required to define 7 macros. 
There were no errors, warnings or information messages. 


| 
CON LCL 
Q.) 01 NOPIC USR CON ABS LCL NOSHR EXE RD — WRT NOVEC BYTE | 
29917.) 2.) PIC USR CON REL LCL SHR EXE RD NOWRT NOVEC LONG | 
MACRO/ENABLE=SUPPRESSION/D1 SABLE= (GLOBAL , TRACEBACK) /L1S=L1S$:BASMATMUL /OBJ=0BJ$:BASMATMUL MSRC$:BASMATMUL/UPDATE=(ENHS: BASMATMUL)+LI_ 
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